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Abstract 

Being able to accurately model and predict the dynamics of dispersed inclusions transported 
by a turbulent flow, remains a challenge with important scientific, environmental and economical 
issues. One critical and difficult point is to correctly describe the dynamics of particles over a 
wide range of sizes and densities. Our measurements show that acceleration statistics of particles 
dispersed in a turbulent flow do exhibit specific size and density effects but that they preserve 
an extremely robust turbulent signature with lognormal fluctuations, regardless of particles size 
and density. This has important consequences in terms of modelling of the turbulent transport of 
dispersed inclusions. 
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Turbulent transport of material inclusions plays an important role in many natural and 
industrial processes. In marine ecosystem for instance, interactions between the turbulent 
sea and the gametes of marine animals is determinant for the efficiency of the reproduction 
and the spreading of the species Industrial stakes concern mixing and combustion, but 
also pollutants dispersion. Being able to accurately model and predict the dynamics of 
dispersed inclusions transported by a turbulent flow, remains a challenge with important 
scientiflc, environmental and economical issues. One critical and difficult point is to correctly 
describe particles' dynamics over a wide range of sizes and densities. The dynamics of 
particles transported by a turbulent flow is known to be affected by size and density effects. If 
particles are neutrally buoyant and small (namely smaller than the smallest turbulent eddies, 
at the energy dissipation scale t] of the flow), they behave as fluid tracers and their dynamics 
reflects fluid particles dynamics. This property is commonly used to characterize single phase 
flows from tracer particles imaging techniques (PIV, LDV, PTV). However, when particles 
density is larger or smaller than the surrounding fluid and/or when particles size become 
comparable to turbulent eddies, their dynamics deviates from that of fluid particles and 
tends to be affected by so called inertial effects. Among them some can be qualitatively 
described in terms of particles interaction with turbulent eddies. For instance, centrifugal 
forces are generally expected to cluster heavier particles in low vorticity regions (and lighter 
in high vorticity regions), an effect known as preferential concentration. However, from 
a quantitative point of view no reliable model has emerged yet to describe and predict 
accurately the statistical properties of particles advected by a turbulent flow. For instance 
we are unable today to predict correctly the dynamics of water droplets in a cloud j^, 
what considerably limits our understanding of rain mechanisms. Actually, even writing 
an appropriate equation of motion for a particle transported by a turbulent flow remains a 
theoretical challenge which has only been approached in some limit cases, assuming generally 
point like particles 

mm 

; but the range of validity of such models for real particles, 
with flnite size and flnite density, as well as the minimal relevant ingredients for a model to 
be pertinent remains unclear. 

In the present study we have measured acceleration statistics of isolated material par- 
ticles transported in a turbulent air flow, varying systematically particles size and density. 
Acceleration is a kinematic quantity of particular interest since it directly reflects the forces 
exerted by the carrier flow on the particles, which is the basic ingredient of any dynamical 
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model for the particles equation of motion. As particles we use soap bubbles ; their density 
can be adjusted from neutrally buoyant to about 70 times heavier than air. As their Weber 
number is extremely small, they are known to not deform and to behave as rigid spheres. 
The seeding density is extremely low (particles are injected individually) so that particles 
can be considered as isolated and do not backreact on the carrier flow. Particles are therefore 
only characterized by two parameters : the ratio of their diameter D to the dissipative 
scale of the flow (0 = D/rj) and the ratio F of their density pp to the carrier fluid (air in 
our case) density p^ir 

(r = Pp/Pair)- Figure [1] summarizes all the particle classes (in the 
(0, r) parameter space) that we have considered in the present study. As far as we know 
this represents the most exhaustive exploration over such a wide range of particles sizes and 
densities ever done in a same experimental configuration. 

Our experiment runs in a large wind tunnel with a measurement section of 0.75 mx 0.75 m 
where the turbulence is generated downstream a grid with a mesh size of 6 cm and reproduces 
almost ideal isotropic turbulence. The results reported here were obtained with a mean 
velocity of the fluid U = 15 m-s"^ and a turbulence level Urms/U ~ 3%. The corresponding 




FIG. 1: Particles classes considered in the present study, described in the (</>, F) phase space. The 
countour lines indicate an estimation of the particles Stokes number St = Tp/ru, where the respone 
time Tp of the particle has been corrected for added mass and for finite Reynolds number effects 



according to Schiller and Nauman 
particle diameter D. 



and T£) is the flow eddy turnover time at the scale of the 



3 




10 



15 



-10 



-sw tfia -1« e . m v» m 



-5 5 



10 



15 



FIG. 2: Normalized component acceleration probability density function (PDF) of material parti- 
cles transported in a turbulent air flow. Dashed line corresponds to the measurement by Ayyaso- 
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mayalula et al. [l3|] for water droplets (with size </> ~ 5 x 10 ^ and T 1000). Dot-dashed line is a 

^n(|x/^/3|)-^2s 
V2s 



fit by the relation 'P{x) 
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the acceleration amplitude 
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associated to a lognormal distribuion of 



15[ (best fit is found for s 0.62, corresponding to a distribution 
flatness J- = |e^*^ ~ 8.4). Inset shows the non normalized component acceleration PDF for F ~ 1 
(outer bundle, green tone) and F '--^ 5 (inner bundle, red tone) particles and various values of (p. 



Reynolds number, based on Taylor microscale, is of the order of Rx = 160. The dissipation 
scale T] = (z/^/e)^/^ is 240 /im (where u is the kinematic viscosity of air and e the turbulent 

'articles 
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energy dissipation rate per unit mass) and the energy injection scale L is 6 cm. \ 
are individually tracked using ID Lagrangian acoustic Doppler velocimetry 

an 

measure the streamwise velocity component Vz of the particles as they are tracked along 
their trajectory. The streamwiase acceleration component is obtained by differentiating 
the velocity using a convolution with a differentiated gaussian kernel 12|. Each particle is 
tracked during approximately 50 ms, which corresponds to several dissipation time scales 

= (z//e)^/^ ~ 3.8ms. For each point in the (</>, F) parameter space we record at least 4000 
tracks at a samplimg rate of 32768 Hz, giving more than 10^ data points per set. 

Figure [2] represents the probability density functions (PDF), of the acceleration compo- 
nent normalized to variance one, (a^/ (o^)^^^), for all studied particles. Remarkably we find 
that, within statistical errorbars in the rare events tails, all the PDFs almost collapse onto a 
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single curve, indicating that normalized acceleration statistics depends significantly neither 
on particle size, nor on particle density. These results extend to the case of heavy particles 
our previous observation for neutrally buoyant particles 15|], where we have shown this PDF 
to be correctly described by a robust lognormal distribution V (yellow dot dashed line on 
figure m see the caption for its analytical expression). They are also consistent with recent 
measurements, in a water von Karman flow, for particles in the range 0.5 < F < 2 

Although for technical reasons, we cannot reduce at present the bubble's diameter below 
~ 10, it is still interesting to compare our results with other recent measurements obtained 
by Z. Warhaft's group at Cornell University for heavy sub-kolmogorov particles (0 <^ 1) 13|. 
Using high speed optical tracking, they have measured acceleration statistics of small water 
droplets (</> ~ 5 x 10~^ ; F ~ 1000 ; St ~ 0.1 ) transported in a turbulent air flow with 
comparable characteristics to ours (both are windtunnel experiments with similar isotropy 
levels), although at a slightly higher Reynolds numbers (i?A ~ 250). We have superimposed 
on figure [2] the normalized acceleration PDF from their measurements (black dashed line). 
We find it to be almost undistinguishable from our measurements with larger particles. 
This suggests that acceleration is a physical quantity with an extremely robust lognormal 
statistical signature, over a very wide range of sizes {D ~ 0.05?7 -^/lO) and densities 
(F ~ 0.5 1000). Once normalized to variance one, all PDFs are well fitted by the same 
distribution V which is only parametrized by its flatness JF (in the present experiments 
J^~8.4). 

A deeper insight into the specific effects associated with particles' size and density can be 
obtained by examining acceleration variance. Since the statistics of acceleration normalized 
to variance one is found essentially independent of particles size and density, an effect of 
those parameters can indeed only be expected to affect the acceleration variance itself (a^). 
The inset in figure [2] represents the true acceleration PDF (not normalized to variance one) 
for neutrally buoyant (F ~ 1) and heavy (F ~ 5) particles, with varying sizes. It appears 
that heavy particles have a much narrower and peaked PDF, indicating a gobal decrease 
of acceleration variance with increasing density. The spreading of the curves for a given 
density indicates that size also influences particles acceleration variance. 

A closer analysis of acceleration variance shows a non-trivial dependence with size and 
density. We report in figure [3] the acceleration variance for all the measurements we have 
performed (in the following, we consider the dimensionless acceleration variance normalized 
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FIG. 3: Acceleration variance for all particles studied in the (0, T) parameter space (dots). The 
surface represents a rough interpolation based on the available measurements. The dot dashed 
lines materialize the interpolation along the best resolved constant c^-lines and F-lines. For the 
F = 1 set, the dot-dashed line coincide with a 0"^/'^ decay for 15 [l^. 

according to Heisenberg-Yaglom's scaling : Ao(0, F) = (a^) e~^/^i/^/^). The mapping of 
the (0, F) parameter space that we have been able to achieve, allows us to infer a rough 
interpolation of the evolution of Ao{(j), F), represented by the surface in figureO Though only 
the coarse tendency of this interpolation is relevant, several important qualitative features 
can still be observed. For small and close to neutral particles Aq ~ 2.8, a value consistent 
with previous measurements [l8| and DNS Q, [2^ for particles in the fluid tracer limit. If we 
consider the effect of an increasing density at a fixed particle size (along (p = const ant- lines), 
Aq is always found to decrease and to saturate to a finite limit (noted a[f (0)) for the largest 
densities. This is better quantified in a 2D projection (figure Hj). The (p ~ 16.5 set of 
measurements (which is our most complete set of density effects at fixed size) shows that 
the evolution of acceleration variance with density exhibits two different regimes (see inset 
in figureH]): (i) for low densities Ao(0, F) = ao(0)r" (ao(0) corresponds to size effects for the 
neutrally buoyant case which we have previously studied in [l^) with a ~ 0.6 (ii) for large 
densities Aq saturates to the finite limit a[^(0) (a[^(16.5) ~ 0.7). The transition between 
these two regimes occurs for a characteristic density ratio F*(0) (F*(16.5) ~ 10). As shown 
on figureH] from the measurements at different values of (p, the existence of these two regimes 
seems to hold for all studied particle sizes but with a size dependent transition density ratio 
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r*(0) and a size dependent saturation value a^(0). If we consider more closely size effects at 
fixed particle density (along F = constant-lines), as seen on figure [3|, the scenario is actually 
rather complex, since depending on the density ratio F, Aq can either decrease or increase 
with particle's size. For neutrally buoyant particles (F = 1) Aq starts to deviate from the 
fluid tracer value for particles larger than about ~ 15 and then decreases monotically 



as We have shown in a previous study 15[ that this is the expected scaling for 



inertial range sized particles when the main forcing simply comes from the spatial pressure 
differences of the unperturbed flow around the particle. As we move to larger density ratios 
we then observe a continuous transition toward a drastically changed size dependence for 
the largest densities (F > F*) where a'^{(j)) experience a sudden increase for sizes around 
17 < < 19. Outside this transition region (i.e. for < 17 and > 19), doesn't 
exhibit significant dependence on (at least in the accessible range of sizes) as also seen 
on figure m where we have a[f (13) ^ a[f (16.5) and a[f (19) ~ a^(21). Inset in figure H] also 
suggests that F*(0) might exhibit a similar transition in the same range of sizes as we have 
F*(13) ~ F*(16.5) ~ 10 and (F*(19), F*(21)) < 5. A size dependence of the exponent a (for 
F < F*(0)) cannot be excluded, though a more detailed mapping of the parameter space in 
the region F < F*(0) is still required to be conclusive. 

Our measurements have important consequences in terms of the development of accurate 
models for the turbulent transport of finite size material particles. The robustness of normal- 
ized acceleration PDF, regardless of particles' size and density, imposes a strong constraint 
on the statistical properties of the forcing terms to be included in the equation of motion. 
For instance, models based on the Maxey & Riley equation 3, J] for inertial point particles 
(characterized mainly by their Stokes number) predict a monotonic decrease of acceleration 
flatness, with a continuous trend of acceleration PDF to gaussianity, as well as a monotonic 
decrease of acceleration variance with increasing particles Stokes number ^, [l3] (which in 
these models is only parametrized by an increasing particle's response time and eventually 
interpreted as an increasing density ratio and/or particle size). These trends are contradic- 
tory with the experimental evidence we report here (acceleration flatness remains size and 
density independent and acceleration variance may increase with particle size in the limit 
of high density ratios). Our results indicate that the high Stokes number extrapolation of 
such point particle models is not adequate to correctly describe the dynamics of finite size 
inertial particles. In other words, inertial range sized particles cannot be correctly modeled 
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by simply considering increasing response times in the point-particle approximation; density 
and finite size effects require a simultaneous specific modeling. 

With this goal, finally, we briefly discuss and propose a simple phenomenology for the 
observed trends of acceleration variance with size and density in the context of the recently 



introduced sweep-stick mechanism 2l|, which offers an interesting frame to interpret 



several of our observations. Physically this mechanism relies on the simple fact that inertial 
particles reside longer in the quietest regions of the flow and tend therefore to stick pref- 

y leads 



23 



2^, 



erentially near zero-acceleration points of the carrier turbulent field (this eventual 
to the so-called preferential concentration effect, responsible for particles clustering 
which in this context is not only attributed to centrifugation by high vorticity regions) along 
which they tend then to be advected. This phenomenology is consistent with the decrease 
of acceleration variance that we observe when particles density increases. Another consis- 
tent point is the fact that velocity statistics of our particles (not shown here) are found 
identical to that of the carrier flow (obtained from classical hotwire anemometry), a feature 
also shown numerically for zero- acceleration points {22 1 (but which is at odds from usual 
predictions based on point-like particles approaches 25|]). Pushed further, the sweep-stick 
phenomenology also offers a consistent frame to possibly explain the sudden increase of ac- 
celeration variance with particle size observed for high density ratios: while heavy and small 
enough particles can indeed be expected to efficiently hide in the quietest regions of the flow, 
when particles become larger than the typical size L* of these regions, the quietening effect 
is damped as particles experience again the influence of active regions in the turbulent field. 
Though the existence of such a typical scale is still controversial, recent measurements of 
preferential concentration suggest that the quiet sticking regions might have a characteristic 



2J] consistent with the 17 < < 19 range for which we observe 



size L* in the range 10 — 20// 
the sudden increase of a[f (0). 

To summarize, using a versatile material particle generator we have been able to explore 
the simultaneous influence of size and density ratio on the turbulent transport of material 
particles. Their acceleration statistics are found to be robustly described by a lognormal 
distribution, where only the variance depends significantly on particle's size and density. We 
have shown that, for heavy particles, finite size effects can be trivially extrapolated neither 
from the heavy point particle case nor from the finite size neutrally buoyant case. However, 
the simultaneous influence of density and size may be consistently interpreted in the context 



8 



3 
2.5' 
2 

< 1.5 
1 

0,5 



10 20 30 40 50 60 70 

r 

FIG. 4: 2D-Projection on the (F — Aq) plane of measurements data points on figure El showing 
the evolution of acceleration variance with density ratio F for different particle sizes (</> ^ 16.5 
set (circles) is the best resolved). The inset shows the same data in a log-log plot. Errorbars are 
mostly due to experimental uncertainties in the determination of particles size and density. 

of sweep-stick mechanisms, what should stimulate further numerical investigations of 
acceleration field in turbulent flows. It is our hope that these measurements will contribute 
bridging experimental observations and developments of accurate models for the turbulent 
transport of particles. 
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